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Статистический анализ систем с цифровыми 
водяными знаками 


В работе рассмотрен статистический подход к построению и оценке качества модели стеганографической 
системы с цифровыми водяными знаками. Описан мегод качественной оценки детектора и декодера ЦВЗ, 
позволяющий гарантировать определенный уровень надежности. Показан путь построения оптимального 
детектора и декодера для систем маркировки, базирующихся на методике расширения спектра сигнала. 
Даны соотношения статистических решающих тестов для детектора при известном вероятностном 
распределении множества возможных изображений и при приведении гистограммы изображений к 


фиксированному виду. 


Введение 


Представление, хранение и распределение информации в цифровом виде давно 
стали повсеместно используемыми, привычными и удобными манипуляциями над 
ней. Вместе с тем легкость восстановления и коррекции данных, возможность 
использования универсальных средств работы над ними и другие преимущества 
цифрового представления информации могут нивелироваться простотой получения и 
распространения копий цифровых объектов, полностью идентичных оригиналу. Это 
влечет за собой необходимость активных исследований проблемы защиты авторских 
прав на цифровую интеллектуальную собственность. 

В связи с недостатком теоретических наработок по проблеме систематизации и 
оценки качества существующего многообразия методов и алгоритмов защиты авторских 
прав в статье показан пример моделирования этой проблемы для схем с аддитив- 
ными цифровыми водяными знаками (ЦВЗ), базирующихся на методике расширения 
спектра сигнала [1], [2]. 

В работе рассматривается статистический подход к построению математической 
модели систем с ЦВЗ, которая в дальнейшем может служить основой для практи- 
ческой разработки эффективных детекторов и декодеров водяного знака. Статья 
посвящена разработке модели детекторов и декодеров, требующих на входе не 
наличия оригинального изображения, а только определенных статистических данных 
о множестве возможных оригинальных изображений (распределение этого мно- 
жества) или о виде гистограмм изображений. 

Такой подход удобен при оценке качества работы системы с ЦВЗ и для 
сравнительного анализа различных методов и алгоритмов маркировки изображений. 
Кроме того он позволяет априори задавать определенный уровень качества путем 
наложения ограничений на допустимый процент ошибок при обнаружении и 
декодировании водяного знака. 
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Общая модель системы маркировки изображений 
и ее характеристики 


Пусть владелец авторского права имеет некоторый секретный ключ К, отоб- 
ражающий его право собственности, и хочет защитить свое цифровое изображение 
Х путем внедрения в него сообщения У, состоящего из элементов некоторого 
конечного дискретного алфавита и однозначно идентифицирующего владельца. 
Рассмотрим данную задачу, опираясь на блок-схему общей модели системы марки- 
ровки изображений (рис. 1). 


Обнаружение 
К т К 


Внедрение ЦВЗ 


Извлечение 


ЦВЗ К 


Рисунок 1 — Блок-схема общей модели системы маркировки изображений 


Согласно схеме оригинальное изображение Х преобразовывается в маркиро- 
ванное изображение У с помощью функции №, на вход которой также подаются ключ 
К исообщение У . Далее изображение У пересылают получателю по открытому каналу 
связи, где оно может подвергнуться естественным искажениям, присутствующим в 
канале, или преднамеренным атакам на ЦВЗ. Причем любые преднамеренные атаки 
в данном случае выполняются без знания К и Х, так как эти элементы публично 
недоступны. Таким образом, изменения У можно промоделировать как канал с 
помехами, вход которого У и выход / связаны условным распределением Руу. 


В процесс проверки авторского права вовлечены две функции: детектор и 
декодер ЦВЗ. Детектор ЦВЗ 4 решает, содержит ли изображение 1, водяной знак, 
сгенерированный при помощи ключа К ‚ т.е. на его вход поступают 7, и К, а на 
выходе содержится булево решение. Если детектор определяет наличие водяного знака 
в изображении, то этим будет доказано авторское право человека, который владеет 
ключом К. В дальнейшем возможно извлечение скрытого сообщения с помощью 
декодера ЦВЗ е, на вход которого подаются 1, и К, а результатом работы является 
оценка скрытого сообщения У . Отметим, что ни у детектора, ни у декодера нет дос- 
тупа к оригинальному изображению Х. 

Качественная оценка детектора ЦВЗ выполняется на основе двух критериев — 
вероятности ложной тревоги Р; и вероятности обнаружения Рь [3]. Р, определяет 
вероятность получения позитивного результата при тестировании 7, на обнаружение 
ЦВЗ, тогда как фактически это изображение не содержит водяного знака, сгенери- 
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рованного с ключом К. Р, - вероятность получения позитивного результата тогда, 
когда изображение действительно содержит ЦВЗ, сгенерированный с ключом К. 
Пусть Н, - это гипотеза «изображение содержит ЦВЗ», а Н, - гипотеза «изображе- 
ние не содержит ЦВЗ». Тогда 


р, =Ри{ (2, К) = Н, Но}, (1) 


ЧЕ 
Р,‚ =Рг{9(7,К)=Н,|Н,\. (2) 
Качественная оценка декодера ЦВЗ - это вероятность ошибки Р‚, определенная 


как вероятность получения неправильной оценки У: 


Р. РУ У}, (3) 


Соотношения (1) — (3) назовем мерой производительности системы маркировки 
изображений. Отметим, что при практической реализации детекторов и декодеров 
особый интерес представляют две следующие условные вероятности ошибок: 


Р.О =РУ = УХ}, (4) 


Р.(К)=РУ #У| К} (5) 


где Р.(Х) - вероятность ошибки, обусловленная для данного оригинального изображе- 
ния Х , Р. (К) - вероятность ошибки, обусловленная для определенного ключа К. Первая 


вероятность показывает, насколько данное изображение подходит для сокрытия в 
нем сообщения, а вторая — среднюю производительность для конкретного владельца 
авторского права. 


Система маркировки, использующая методику 
расширения спектра 


Многие из существующих методов и алгоритмов создания ЦВЗ так или иначе 
эксплуатируют методику расширения спектра или, другими словами, широкополос- 
ную модуляцию. Методика широкополосной модуляции позаимствована из теории 
связи [1], [2]. Дело в том, что сокрытие данных можно рассматривать как ком- 
муникационную проблему, в которой оригинальное изображение играет роль шума 
канала, а атакующие могут попытаться разрушить передачу информации. 

Пусть необходимо скрыть М информационных бит в изображении Х = (х, КЕ ) . 
где Г.- количество элементов в нем. Для этого сначала с помощью генератора 
псевдослучайных чисел (ПСЧ), инициализируемого с помощью ключа К, образу- 
ется псевдослучайная последовательность $= (61,5....31). Далее для того, чтобы 


гарантировать невидимость водяного знака, 3 поэлементно перемножается с перцепци- 
онной маской б [4]. Чтобы получить перцепционную маску, используют психови- 
зуальную модель, которая строится на основе учета свойств системы человеческого 
зрения (СЧЗ). Перечислим некоторые из этих свойств: 

— яркостная чувствительность (СЧЗ более чувствительна к шуму на участках с ма- 
лой яркостью); 
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— частотная чувствительность (глаз более восприимчив к низкочастотному, чем к 
высокочастотному шуму); 
— эффект маскирования (аддитивный шум гораздо заметнее на гладких участках 
изображения, нежели на высокочастотных); 
— чувствительность к контрасту (глаз наиболее чувствителен к изменениям высоко- 
контрастных участков); 
— чувствительность к размеру (большие детали заметнее); 
— чувствительность к цвету (некоторые цвета, например, красный заметнее других); 
— чувствительность к местоположению (человек обращает больше внимания на 
центральную часть изображения) и т.д. 

Построенная на основе вышеперечисленных свойств маска б учитывает, как 
изменение отдельных элементов изображения влияет на его полное перцепционное 
искажение. То есть с помощью 6 добиваются максимальной визуальной незамет- 


ности ЦВЗ. Набор индексов НВ. разделяют на подмножества ©, 


1? 


так что 
УТ ] ООО, =@. В дальнейшем в каждый о, будет внедрен ]-ый бит 


сообщения. Для обеспечения дополнительной стойкости к атакам такое деление 
может быть зависящим от ключа. Водяной знак генерируется согласно правилу 
\,; =У,6:5;, МЕО,, (6) 


где ] Е {,2. з №}, а коэффициент у, кодирует ]-ый бит скрываемого сообщения. 


Перцепционный Широкополосная 
анализ модуляция 


Рисунок 2 — Схема маркировки изображения, использующая методику 
расширения спектра 


В матричном виде 
\=А(К,ХУ, (7) 
ЧЕ В 
где У=(у,,Уу....ух) — скрытое сообщение, а А(К,Х) -— матрица размером [Ех М, 
элементы которой удовлетворяют условию: 
6:5; 1Е О, 
а, = . у (8) 
0, 1#0, 
Столбцы матрицы А =(а,,а,...ах) назовем модуляционными импульсами, так как 
цифровой водяной знак может быть выражен как их линейная комбинация, что мож- 
но сравнить со схемами мультиимпульсной амплитудной модуляции [5], [6]. 
Итоговое маркированное изображение формируется как сумма оригинального 
изображения и водяного знака У =Х+\У\ (рис. 2). 
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Таким образом, в процессе генерирования ЦВЗ вектор скрываемого сообще- 
ния, определенный в М-мерном пространстве, случайным образом отображается на 
пространство большей размерности (Г, >> М), что в теории связи и называют 


методикой расширения спектра (Зргеа4 Зресгат). 


Получение изображения с заданной гистограммой 
(плотностью распределения) 


Источником получения полезной статистики об изображении является гистограм- 
ма. Перед тем, как выполнять маркировку изображения согласно схеме, описанной в 
предыдущем разделе, мы предлагаем привести его гистограмму к заданному, заранее 
определенному виду. Более того, если известно множество всех возможных ориги- 
нальных изображений, мы предлагаем выполнить приведение гистограммы для 
каждого изображения из этого множества. 


Пусть дано некоторое оригинальное необработанное изображение С = (2 а ). 
Яркость пикселей изображения 1(С) является дискретной случайной величиной с 
функцией распределения Р(С)= {Р, (С)}, где Р.(С)- оценка вероятности появления 
пикселя со значением яркости [;,, 1е {0...М _ 1, М - количество уровней яркости. 
Функцию Р.(С) называют также нормализованной гистограммой цифрового изобра- 
жения С: Р/ (Ц) = ее : УР.(6) =1, где г.(С) — количество пикселей с яркостью 

1=1 
Г., Е — общее количество пикселей в изображении. 

Изображение С необходимо преобразовать таким образом, чтобы его яркость 

имела априори заданное распределение Р(Х)= {р СО ] Е {0...М = Будем счи- 


тать, что уровни яркости упорядочены по возрастанию. 
Алгоритм отображения Р(С) в Р(Х) имеет вид: 


Шаг 1. 
1) если Р(@) =Р‚(Х) , то ничего не делаем; 


2) если Р(@)>Р,(Х), то количество пикселей в исходном изображении г, (С) 
нужно уменьшить до величины г(Х) путем увеличения яркости 
|“ (ХХ -к, () пикселей; 

3) если Р.(@)<Р.(Х), то количество пикселей в исходном изображении 
г (С) нужно увеличить до величины г,(Х) путем уменьшения яркости бо- 
лее ярких | (Х) (С) пикселей в С до |; 

4) пересчитываем вероятности Р;(@) после сделанных преобразований. 

Шаг 2. Имеем Р‚(@)-+Р! (С) =Р,(Х)+Р,(@) =Ре(Г), 

РОО -Р:(Х) =Рх (1). 

1) если Р/(@) =Р!(Х) ‚ то ничего не делаем; 

2) если Р(9)>Р!(Х), то количество пикселей в текущем изображении с яр- 
костью [ нужно уменьшить путем увеличения яркости | (Х-п (С) пик- 


селей; 
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3) если Р,(@) <Р/(Х) ‚ то в изображении @ выбирается | (Х)- (С) пикселей 


с яркостью, большей нежели [| и их яркость уменьшается до [; 
4) пересчитываем вероятности Р; (Ц). 
Третий и все последующие шаги выполняются аналогично второму шагу. 


Статистическая задача о выборе из двух гипотез 


Рассмотрим функционирование детектора ЦВЗ как решение статистической 
задачи проверки основной гипотезы Н, -— «изображение не содержит ЦВЗ», при 
условии, что для нее имеется всего одна альтернатива, гипотеза Н,-— «изображение 


содержит ЦВЗ». Цель раздела состоит в описании критерия выбора из двух гипотез, 
что в дальнейшем позволит конструировать оптимальные тесты на обнаружение 
водяного знака. 

Всякое правило принятия решения характеризуется вероятностью принять ту 
или иную гипотезу в зависимости от наблюдаемого значения случайной величины 
7. Пусть п(7.) есть «критическая» вероятность — вероятность отвергнуть основную 


гипотезу, если наблюдаемое изображение есть ,, т.е. 
п(7) =Р(Н, | 2) (9) 


(вероятность принять гипотезу Н, о наличии ЦВЗ в изображении при условии, что 


наблюдалось 7.). Качество правила определяется вероятностями принятия и отвер- 
жения каждой из гипотез в зависимости от того, какая из гипотез верна. Его 
характеризуют вероятностями ошибок. 
Ошибка первого рода — отвергнуть истинную основную гипотезу Н,, она 
обычно обозначается © 
о =Р(Н, |Н,). (10) 


Ошибка второго рода — принять основную гипотезу, если верна альтернатива, 
ее вероятность обозначается В: 


В=Р(Н,|Н)). (0 


© — уровень значимости критерия, 1-В — его мощность [7]. Согласно введенным в 
первом разделе обозначениям © =Р, (уровень значимости соответствует вероят- 
ности ложной тревоги для детектора ЦВЗ), 1-В=Рь (мощность — это вероятность 


обнаружения ЦВЗ детектором). Уровень значимости выбирается заранее. В нашем 
случае при выборе нужно исходить из коммерческой целесообразности допусти- 
мости ложной тревоги для системы с ЦВЗ, например, уровень значимости 1 %, при 
этом а = 0,01. Мощность критерия максимизируют, т.е. вероятность ошибки второго 


рода стараются сделать минимальной. Если выбрана критическая вероятность, то 
а = [92 (47), В = Га —п(7))Е (47), (12) 
где Е (7) — распределение величины при гипотезе Н; (1 = 0,1). 


Наиболее «мощный» критерий с заданным уровнем значимости — критерий 
Неймана-Пирсона. Он определен для случая, когда мера Е абсолютно непрерывна 
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ш 


относительно меры Е: для всех борелевских множеств СЕ К 


в(©= [28 2). (13) 


Рассмотрим два важных случая, когда равенство (13) выполнено. 

1) Пусть ЕЁ и Е имеют плотности относительно лебеговой меры в К", равные 
Е (2) и Е,(2) соответственно, причем Ё, > 0. Тогда равенство (13) выполняется, 
если 


Ё 
Е(7) = Е те | (14) 


2) Пусть обе меры Е(47,) дискретны, т.е. можно указать такую последо- 
вательность 7, ЕВ", что УЕ, (17, } = УЕ, ({7, }) =|иЕ ({7.}) > 0 для всех 7,. Тогда 


ое из 


(значения #(7,) при 7 * 7, роли не играют). 


Пусть /, — такое множество, что #(7)>1{ при ХЕХ, #(7) <{ при ДЕК" -Т.. 
При Е (7,) = © критерий является наиболее «мощным» при уровне значимости © для 
непрерывного случая. Кроме того всегда можно разделить В" на два множества 
Ри В" -БО,„ такие, что #(7)> 1, на), и Ё(7)< Е, на В" -ОБ, и Е(Ф,„)=а. 

Для дискретного случая К" делится на три множества О „, где #(7)> 1, Г, 


где Е(7) =а, и В" -Р, Г. Если л(7.) для критерия выбирается как 


17 ЕО 
п(7,) = р, ЕГ, ‚ 
0,7 ЕВ" -О„, Г, 


о=рЕ (Г. )+Е (О, ), 


то критерий также является наиболее «мощным». 


Построение оптимальных структур детектора 
и декодера ЦВЗ как решение статистической задачи 


Пусть имеется изображение 1, и ключ К ‚ а также предполагается, что марки- 
рованное изображение не переносило ни неумышленных искажений, ни атак, т.е. 
7=У. Цель теста на обнаружение ЦВЗ состоит в том, чтобы решить, содержит ли 
изображение 7, водяной знак, сгенерированный владельцем авторского права, кото- 
рый обладает ключом К. Этот тест можно сформулировать как бинарный тест двух 
гипотез: 


Н,: (=Х, +А(К,Х, У, 
Н.: 7=Х,, (16) 
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где Х и Х - какие-то изображения. Так как декодировать скрытое сообщение в 
1 0 


данном случае не обязательно, элемент У рассматривается как случайный вектор с 
функцией распределения, равной вероятности распределения сообщения. 

Чтобы убрать зависимость А от неизвестного детектору оригинального изобра- 
жения, А(К,Х) разумно аппроксимировать заменой на А(К, 7), так как ЦВЗ при- 


вносит в изображение небольшие искажения, которые ожидаемо повлекут за собой 
незначительные изменения в перцепционной маске 6. Тогда тест (16) можно 
аппроксимировать следующим бинарным тестом: 


Н,: 7=Х, +А(К, ИУ, 
Н.: Й=Х.. (17) 


Пусть ЗЕ НН} — решение, полученное при тесте на обнаружение ЦВЗ (фун- 
кция 4 на рис. 1). Чтобы гарантировать надежность системы маркировки изобра- 
жений, вероятность ложной тревоги, определяемая как Р; = Рг{$ =Н, |Но}, должна 
быть ниже некоторого максимального значения (порога). При разработке детектора 
преследуется цель максимизировать вероятность Рь = Рг $ =Н, |1}, которая соот- 
ветствует максимально допустимому Р,. 

В тесте на декодирование ЦВЗ предполагается, что 7, действительно содержит 
водяной знак, принадлежащий обладателю авторского права, владеющему секрет- 
ным ключом. Цель декодера состоит в том, чтобы получить оценку \ вектора сооб- 
щения У таким образом, чтобы минимизировать вероятность ошибки. Для декодера, 
как и для детектора, матрицу А(К,Х) можно аппроксимировать матрицей А(К, 7). 

Принимая ключ К фиксированным, можно сказать, что в статистических ре- 
шающих тестах есть два случайных вектора — Х и У. Если известно распределение 
Р.С), оптимальные решающие тесты образуются согласно правилу Неймана- 
Пирсона, рассмотренному в предыдущем разделе. Например, в тесте на обнаружение 
ЦВЗ оптимальный детектор, который максимизирует Рь для любого значения Р,, 
задается тестом: 


НЬЮ) _ пу. -АЬВУ) > 


вен 5 Е, (2) Г. е) 


0 


где п- порог, величина которого зависит от @а и распределения Е (7) согласно (12). 


1 
Если принять сообщения равновероятными, то УУ р(У) = м? где М - количество 


элементов алфавита, из которого они формируются. 
Оптимальный декодер ЦВЗ, который минимизирует вероятность ошибки, 
обусловленной для ключа К ‚ задается структурой максимальной вероятности: 


У =аг8 тах шЁ, (7 | У) = агё тах ш Е; (7 АУ), (19) 


где приставка аго означает, что ищется не максимальное значение, а номер макси- 
мального аргумента функции тах. 
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Выражения (18) и (19) задают функции детектора и декодера, представленных 
на рис. 1, 4(7,К) и е(7,К.) соответственно. 


Если распределение Ёх(Х) не известно, но известно, что для множества всех 
возможных оригинальных изображений их гистограммы имеют одинаковый вид — 
Р(Х) = {р (х)} ] в {0л...М 1}, то критерий Неймана-Пирсона следует использовать 
для каждого уровня яркости. Если хотя бы для одного уровня яркости принимается 
решение в пользу гипотезы Н,, то детектор принимает решение о наличии ЦВЗ. 
Оптимальные решающие тесты для каждого уровня яркости ] получаются из условия 


НЕЮ В 
Р.(2|Н.) Р.Х < 


(20) 


где ^. — порог, Р.(Х) и Р‚ (7) — вероятности появления }-ого уровня яркости в ори- 


гинальном и наблюдаемом изображениях. 


Влияние искажений и атак 


В процессе распространения маркированное изображение может быть случайно 
или специально искажено, например, аддитивным шумом, фильтрацией или геометри- 
ческими преобразованиями. В целом множество возможных искажений ограничивается 
искажениями, которые не приводят к чрезмерной деградации изображения, делающей 
его непригодным к использованию. Искажения и атаки ухудшают производительность 
решающих тестов. Так как все возможные деформации изображения смоделировать 
трудно, детектор ЦВЗ следует разрабатывать так, чтобы максимизировать Рь для самого 


плохого случая, связанного с атакой самого плохого случая. Аналогично декодер следует 
разрабатывать так, чтобы минимизировать вероятность ошибки Р, для самого плохого 


случая. 

Самая плохая атака для рассмотренной системы маркировки — атака при помощи 
геометрических преобразований, таких как сдвиг, масштабирование, вращение, обрезка 
изображения и пр. Чувствительность детектора ЦВЗ к этому виду манипуляций про- 
истекает из «белой» природы псевдослучайной последовательности $ = (5,5 те ь}} 


что приводит к десинхронизации модуляционных импульсов, фактически представленных 
в изображении и сгенерированных во время решающего теста. Перспективным под- 
ходом для достижения надежности против геометрических искажений является ком- 
пенсация геометрических преобразований перед извлечением водяного знака либо же 
внедрение ЦВЗ в инвариантную к геометрическим преобразованиям область [8], [9]. 


Заключение 


В дальнейшей работе планируется продолжить статистический анализ систем 
маркировки изображений. В частности, для систем, использующих методику расшире- 
ния спектра, рассмотреть случай построения оптимального декодера при неизвест- 
ном распределении оригинальных изображений, но с априори известной единой 
гистограммой. А также случай, когда гистограмма оригинального изображения не 
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известна, но известна другая статистическая информация, такая как, например, 
дисперсия и матрица ковариации [10]. Кроме того планируется рассмотреть при- 
менение кодов коррекции ошибок в процессе генерации ЦВЗ, что позволит повысить 
надежность системы маркировки. 
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В.К. Задрака, Н.В. Кошюна, Л.Л. Яютенко 

Статистичний анал! систем з цифровими водяними знаками 

У робот! розглянуто статистичний шдхд до побудови та ошнки якост! модел! стеганографчно! 
системи з цифровими водяними знаками. Описано метод як!сно! ощнки детектора та декодера ЦВЗ, 
який дозволяе гарантувати певний рйвень надйност!. Показано шлях побудови оптимальних 
детектора та декодера для систем маркування, як! базуються на методиц! розширення спектра 
сигналу. Подан! сшвв1дношення статистичних вирипуючих теств для детектора при в1домому 
1мовфрненому розпод1л1 множини можливих зображень та при наведенн! г1стограми зображень до 
факсованого виду. 


И.К. Хафгака, М.Г. КозйКюпа, Г..1.. Мкиепко 

'Тве УЗай$са! Апа[уз$1$ о? пе УУжегтагК Зу$фетб$ 

ТБе рарег 15 деуое4 уу арргоасБ ю Фе то4деПп® ап4 Фе диаШу гае оЁ Фе %егапостарВу зузет то4е| 
у Фе @оца| ууайеглагК$. Те диаШу гаёе тео4 оЁ Фе даесюг ап4 десо4дег оЁ фе Иа] ууаёегтатК 1$ 
ейпе4 фах 1$ регпуз1уе зрес1Не4 опагатцее4 гобизтез$ [еуе]. Тре то4еПпе ра оЁ Фе орйта| д&есюг 
ап4 десодег Гог Фе Назоте зузетз 15 геаПхе4 фай \уеге Коип4е4 оп \е $1епа| зргеа4 зресблит тефоч. 
ТБе айзИса| Чес1Чте {е5ё ге!айоп$ Юг фе даесюг аге ргезеще \Теп ргобабИиу ФзелБаноп оЁ Фе 
розз1Ые ппазе 5е{ уаз Кпо\т ап4 \Веп фе ппазе В15юстапаз \уаге Яхеа азШоп. 
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